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GRAVITATIONAL RADIATION IN FLAT SPACETIME RELATIVITY 


J«unes Clark Keith* 

Research Institute) University of Detroit 


ABSTRACT 

The possibility of a 1/c 3 order gravitational radiation is consid¬ 
ered within the framework of the flat spacetime relativistic gravitational 
theory of Gi D» BlrkholT, The method stresses t.ne complementarity between 
results derived by considering direct interparticle actions and the inten¬ 
sity of the far gravitational radiation field. It is proved that centrif¬ 
ugal as well as gravitational forces may dampen the accelerated motion of 
a system of masses. The results predict that a "freely" rotating mass 
system will dissipate energy at a rate directly proportional to its total 
energy and to the factor -k(aw/c) 3 w 1/sec, w being the angular velocity, 
'a' the radius of rotation, c the velocity of light in vacuo, and k a di¬ 
mensionless ratio of the system's moments. 

lor a small rotor magnetically suspended and spun in ultrahigh vac¬ 
uum, the gravitational power loss can be detected directly by measuring 
the angular deceleration. An experimental decision should help shed light 
On the validity of the more generally accepted i/c 6 order radiation pre¬ 
dicted by bther means. 


Research supported in part Ly the U. S. Army Transportation Research 
Command, tort Eustis, Virginia, under Contract No. DA-44-177-TC-f>19i 





Purposei 


b'sing the typical flat spacetine relativistic gravitational theory of (i. D. 
Birkhoff,* we shall apply the theory to prediction of possible gravitational effects 
of an order higher than 1/e . Several effects of interest appear in connection with 
the associated gravitational field and radiation from accelerating masses. 

One result of special interest is the prediction of a new, alternate type of 
•pure' gravitational radiation which lies directly between the 1/c order gravita¬ 
tional aberration effect suggested by Laplace*“ and the 1/c order gravitational 
quadruprtle radiations predicted by Einstein and others.® The new predictions appear 
susceptible to experimental verification. This is interesting in that na gravi¬ 
tational phenomena of an order higher than 1/c have ever been observed. In the 
last section an experlmentum crucis Is proposed to test the validity of the resulting 
predictions. 


Notations needed. 


= gjj - (1,-1,-1,-1) for i - j 6 tj - (1,1,1,1) for i = J 

=Ofori/j ^ 0 for i / J 

i,.1 ,k»m»n= 1,2,3,4 (spacetime components). a = 2,3,4 (space components only). 

□ 2 = fi • □ = g iJ sVaxW wher» the standard summation convention nolds on 
upper-lower repeated indices. V is the 3-vector del operator. 

Y the gravitational constant, Y ~ 1/(1-[v/c] 2 )“ the Lorentz factor, 

h^ gravitational potential tensor. gravitational field tensor. 

T^j, t^(r,ct) stress-energy-momentum tensors of mass and gravitational field, 
u^ 4-velocity component. v 3-vector velocity. 

5 3-vector acceleration, a = a*a 


J 

a) angtili 


4-force component. F 3-vector force. 

elocity (radians./sec). c = 3 x 10 cm/sec speed of light in vacuo. 


Inratiohalized e.g.s. wlien needed, 



I. BAbIS AND SOME INTERESTING CONCLUSIONS Of THE BIRKHOFF THEORY 


I .A Postulates of the relativistic gravitational theory of C* D. Birkhoff* 

The formal basis of the Birkhoff gravitational theory is embodied in the fol¬ 
lowing postulates expounded by Birkhoff, Graef, and Barajas!^ 

fl.l The frame of reference of gravitational phenomena is the flat 
spacetime of Loretitz-Minkowski* 

Bii? The formal principle of equivalence is assumed* 0 

BsfJ The gravitational field is completely characterized in spacetime by a doubly 
Covariaht symmetrical potential tensor hjj. The gravitational field components 
B^jk are linear combinations of the first partial derivatives of hjj t 

ti.4 The h A | of a mass point at rest with respect to an inertial system is equal 
in Spacetime to the product of the Newtonian potential li 0 = yM/r times the 
doubly covariant Kronecker delta. 

B»8 Oravitatiohal perturbations propagate with the velocity of light in vacuo 
ih an ihertial reference system. 

B.6 The h t j of a mass point in arbitrary motion is equal to the set of all the 
instantaneous gravitational fields generated by this mass point in all its posi¬ 
tions in physical space* 

Bit The potential tensor of a gravitational field cadfcerf by many mass points is 
equal to the slim of the gravitational potential tensors of each mass point* 

I should like to draw special attention to Birkhoff*s crucial concept of the 
limiting velocity of ** disturbance" in the perfect fluid of his theory. 0 In light 
of the aklomatie nature of this concept, B.5 might perhaps be worded more suggest¬ 
ively! 

B«ft Inertial and gravitational disturbances, forces, and field potentials prop¬ 
agate at fundamental velocity c of magnitude equal to the speed of light in vacuo* 
independent of the mass or energy density or velocity of the frame of reference* 

Some very pointed arguments of Moshinski^ concerning the complementary pure 
gravitational field interactions can be summarized by adding two postulates to the 
above system* These give a concise description of the field aspects of*Birkhoff 
Theory. 
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H.8 s bach component of the potential tensor hjj in rectangular orthogonal coor¬ 


dinates satisfies the inhomogeneous wave equation: 


□ h i 


! 4nyp 


"1J -A'y 

Tab field iriteraction energy density is given by the Lagrangian function,' 

*L = -1/Rny g 1J dh^/ax 1 'dh w /k J + 4h mn P t rrin ® 
e p r i» the equivalent. mass density of 4he external field. 


Bil - B.7 embody the direct interparticle gravitational action description in 
much the same way the "hew" elcoti«odynaniic« of Moon ancl Spencer 9 gleaned from the 
iBtii Century putties describes the direct, interparticle action in electromagnetic 
theory. Postulates Bt8 and B.0 manifest the field aspects of the gravitational theory 
proposed by Birkhoffi Both field and particle descriptions will be of interest* 

i'he above postulates are consistent with ttirkhoff's initial formulation (1842) 
and must be assumed to generate a completely self-contained flat spacetime relativ¬ 
istic theory of gravitational interaction. The theory as it stands is i strong se¬ 
quel to the first attempts by Nordstrom (1912) and Whitehead (1922) to develop a rel¬ 
ativistic gravitational theory within the uniform flat spacetime of Lorentz-Minkowski, 


t*B Applications ol 1 the theory . 

five gravitational theory of Birkhoff may bo summed up by the equations:* 


c 2 d 2 x^/ds 2 “ f®/m = (3h^/6x k - Oh^dx 1 )^^’ 


3 2 b mrf /3x*3*^ = 4rtyp^ or in free space, 0* 

(i,j,m,ft run from 1 to 4) 


(i) 


The "true" mass m is assumed ««ns,tabt along the world tube n b"; c ,2 J 2 x^/ds 2 = f^/m 
Is the 4-acceleration component ef a particle acted on by the gravitational force, 
and measures the curvature of the world tube which that (extended) particle pervades. 
u-J are components of the world velocity of m and constitute the unit tangent vector 
to the world tube, g*J is the metric tensor for a flat spacetime, hjj are sym¬ 
metric components of the gravitational tensor potential. 

Ih a free stationary space surrounding and proper to a point mass M, the Poisson 
equation of (1) becomes time-independent and yields a proper static solution for li^ji^ 
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^(point mass) 


where is the Kronecker delta, * 1 when i*J and 0 when i^J. This is also Birk- 
hoff's solution for the static, spherically symmetric central mass where integration 
eitends over the volume of "true" mass of the Birkhoff perfect fluid of density 
f>ij • point-mass and spherically symmetric solutions of h^ will be 

considered in this article unless further integration is indicated. 

Graef has shown*® that when the mass M moves with arbitrary velocity, the rel¬ 
ativistic potential encountered in an inertial frame of a remote observer is ob¬ 
tained by ordinary transformation laws for a second order tensor in flat spacetime: 


h ij * (2u i u j " «ij) h o < 8 ) 

Mere U^, Uj are covariant components of relative velocity between M of h Q and the 
rest frame of an inertial observer, h 0 being the Newtonian potential -yM/r. Graef 
also shows how more refined considerations will require all quantities to be taken 
at their retarded times with r in the denominator of h 0 replaced by [r - (?* vj/ej 
in usual manneri 10 

Substitution of the static h of (2) into the force equations (1) gives the 

ij 

equations df motion for a small point mass m moving in the field of M with relative 
velocity components u^* u k * Birkhoff integrated the force equations and derived 
the correct cruciai test predictions from his new theory in i042. The theory so 
far has predicted accurately all observed gravitational phenomena. The correct 
factor of 2 in the gravitational bending of light appears as a direct and non- 
arbitrar^ consequence of the Theory as opposed to such implications as drawn by 
W. Thirring*** In addition Graef has predicted the apsidal advance in the general 
problem Of two bodies without excessive effort,'* and Romero Juarez has done like¬ 
wise With a restricted three-body problem of great interest.* 2 
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* • # 

In 1945, Alba realized that were- a central mass M to rotate about Its axis, its 
gravitational field would be modified somewhat as suggested by the relativistic de¬ 
pendence of the potential tensor In (5). 13 It was not until 1952 that Alba applied 

Ills result to the calculation of non-planar forces on the orbit of a satellite,^ 

© 

thus generalizing Hirkhoff's equations for the relativistic motion of a point mass 
In a static field. 

AH/a*s arguments seem satisfactory from a phenomenological point of view and 
consistent within the framework of the liirkhoff Theory. The axial rotation of an 
extended mass M perturbs its static field because of relation (.T) derived by Uraef, 
The static field of M in its proper frame must first be transformed into a frame ro¬ 
tating with minus the angular frequency of M. Thereby a new field of M is set up 
in all inertial spaces, in particular (to u good approximation) in that space in 
which a satellite orbits about M. Once this new field of M has been found by (9), 
the remaining relativistic corrections depend only on the velocity of the orbiting 
mass m with respect to the center of mass of M. 

Spin acceleration thus appears to perturb the static central field in a singular 
and Irremovable manner. Since Graef has shown that appears diagonalized only 
in its proper rest frame, it can be inferred that, absolute motion or acceleration 

» 

of a source M will in general always affect its external field in the non-accelerated 
"free" space posed by the universe, as will be shown more rigorously in a later 
section. 

The external gravitational fields of accelerating masses are therefore distin¬ 
guishable in the liirkhoff Theory in an absolute sense. The new external field 
potential of ah accelerating mass has the character of a non-diagonalized tensor 
of second order. The complexity of the total external field is thus enhanced in a 
definitive way. 

A natural extension of Alba's argument treats the effect of acceleration of 
the orbital mass m. Let v'^ be components of 4-veiocity of M relative to the center 
of mass of an m-M system, and u*^ be the components of velocity of m also relative 
to the system's center of mass. Let us compare the results of substituting tne 
potentials gi\en t>> (:’) for m into the force equations fl). The acceleration 





* 0 . 

experienced by U as a result or the force exerted on it by the field of m follows 
immediately: 

cVxj/ds* - < 2 u 1 m j vJ “ v 4 )v k 'dh 0 /bx k - <2UjV J u k v k - 1) Ghj/dx 1 . 

* 9fc o %\ k <9u 1 *j/0x ,C - dUjU^/Bx 1 ) (4) 

The last tem drops out umly it the non-diagohal couponents of the tensor potential 
vanish• The acorn1eration expwrienred |>y to as a result of the force exerted on it 
by M obtains by permuting the role* of the v4, u^, and n 0 in (4) . 

To illustrate the physical signIfleanee of these equations, consider the 
ideal situation where » rotates in a circle about M so that h Q and y m 1/ll - (v/c) ] 
are time-independent cvnstants of the motion and the ratio ■/>«! • Under these 
conditions the 3-vector acceleration »1 M is found from equation (4) to be: 

Fg/M = - l)^l 0 - 2 y 2 (2y* -l)h 0 a + /c* (5) 

where h 0 • “Y»/r and «* is m*s acceleration in tlio positive r direction. Hie last 
term of ifi) shows the gravitational Torce acts also on the equivalent mass of the 
m-M lield energy* The lactor ? suggests the gravitational mass is twice the equiv¬ 
alent field energy mass* which appears to be localized about m. When m accelerates, 
tiie equivalent mass of its field interaction energy also accelerates. The first 
term oh i*ie right of. (0) is ttie main relativistic correction term of the BirHioff 
Theory. If the ratls of m /M is not vanishingly small, additional residual terms 
will appear in the description of the motion because of the velocity of M relative 
to the center of mass. 

I.C. Question o f t heoret ic al are! experimental c omparisons 
between tlie grr.v 1 tat tonal theories of Einstein and dirkhofl. 

It would be interesting to fi.id reasoiable in-range experimental tests lor de¬ 
ciding definitely In favor of tne curved or flat spacetime versions of gravitational 
theories. One wonders from the standpoint of the models involved how meny residual 
levels of fine effects the r.insteln and birklioff theories in fact describe. 
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Verification of 1/c 2 order "crucial tests" such as advance of Mercury's peri-^ 
nelion, bending of light around the sun, and shift of the emission line of Mdssbauer 

Y-rays in the Earth's gravitational and rotary centrifugal fields, brilliantly con- ® 

« 

firm the Birklioff theoretical predictions as well as those of Einstein* At the same 
time, the M&ssbauer experiments check the locally indistinguishable nature of the ef¬ 
fects of gravitational and centrifugal fields on electromagnetic waves as commonly 
assumed in. both theories. Such close theoretical correspondence makes one skeptical 
that decisive tests to decide in favor of either theory will by found in the area 
of 1/c* effects, The theories down to this level appear equally valid representations 
of nature. 

This' situation suggests the desirability of examining residual gravitational 
effects of order higher than 1/c*. Do .accelerating systems radiate away energy 
through their gravitational field and run down? • 

There are yet no experiments which determine the existence, magnitude, or velo¬ 
city of propagation of "pure" gravitational radiation (or even induction) fields. 

The "crucial tests" above indicate very indirectly that the characteristic velocity 
of gravity may well equal the velocity of light in vaewo*^ Accnrding to Laplace, if 
gravity propagates at some fundamental velocity, aberration of the central farces ex¬ 
perienced by a satellite should cause 1/c order tangential forces to afibct continue 
ously the satellite's secular erfUation.^ Paradoxically, If c is the speed of light 
in vacuo, such a first order effect would have been observed, astronomically long ag#. 

On the other hand, the l/c e order Einstein gravitational radiation reaction 
(which, Eddington considers a residual, fifth order Laplace effect) remains beyond 
the tdllge of practical measurement 

The direction fakei* by the previous considerations o£-. Birkhoff, Uraef, Moshinski, 

. ’ I 

and Alba have led me to apply the Birkhoff Theory for the first time to effects of an ' 

ofder higher than 1/c*. In the absence of experimental results, such residual grav¬ 
itational effects are only possible descriptions of the physical systems involved. 

This effort Is motivated in part by the fact that the ensuing predictions appear barely 
to be within the range of experimental observation. If this were Indeed the case, the . 
results offer a definitive experimental test for the respective models deriving from 
Einstein and Birklioff theories for the cases of gravitational radiation from acceler¬ 
ating mass systems. 
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Hi liAiilA’lION RLALIION ON MSS AtCLimTING IN GKAVITAflONAL 1IELDS 

i 

II.A Time component of gravitational 4-force on a mass m. 

!o study in a simple way the gravitational radiation process, it would be well 

to formulate e measure of tne total time rate of change of energy in a mass system. 

® o 

With only gravitational forces acting, the time component of inertial force according 

to the Birkhoff theory must equal the time component of gravitational force, i.e. 


= 

Yd/dt(ymc^ = cf* 

(6) 

By imposing Birkhoff* condition 

of orthogonality of the 

world velocity and world 

force, the corresponding relation 

u'^f| = 0 gives that 



lu’fj + u*rf ♦ u 4 fj) 

(7) 

or 



tf* 

= -ef* = -v • 1 

f?a> 


R R . 


Where x is the nwn-relativistic 0-vecior velocity and F^ is the relativistic Birkhoff 
ft-vector force given by summing the spatial <3,.",4) components of the gravitational 
4-force ddfined by II). Combining (f»l and <7a) gives t.lte net time rate of change 
©f energy of the siass m in an external gravitational f ield? 

yil/dtiymc^ = v • = dE/dt (R* 

average of huts cuwmnity over a long t#«e *># |s»tegtal number of revolutions J 1 
is defined by the time-averaging operation: ® 

dE/dT 1/1 ( T v * F^dt, 9 (9) 

' o 5 

With a long term or integral cycle gravitational radiation loss, ttie time component 
of gravitational force would not average to zero and would affect m’s inertial 
energy. Since the gravitational induction field is usually considered the conser¬ 
vative part of the total field, its fluctuations average to zero in this process. 







Since the Blrkhoff Theory is intended to describe completely all gravitational in¬ 
teractions, and the conservative induction field components vanish,* we assume what 
is left from equation (9) contains Information about the gravitational radiation or 
radiation-reaction undergone by an accelerating system of masses• The time-averag¬ 
ing operation on dE/dt in.(8) is implicitly carried out in our following results. 

li.B Subsidiary field conditions. 

Some refinement is necessary to account for the finite velocity of propagation 
of fields. Such subsidiary conditions derive naturally from the postulates of the 
first section. Graef derives retardation rnndltions for the gravitational potential 
tensor using postulates B.0-B.8. 10 Except that they apply to a second order tensor, 
these conditions are tbe same as those of the standard Lien*rd-Wiechert retarded 
potentials of electromagnetic theory so familiar In the literature* Consistent 
results can be obtained Usually by expansion of the fields and velocities measured st 
tbe origin of the observer in s Taylor Series about the point t 0 - r { /c, where rj is 
the distance to each source-mass element dmj. Ihiyslcally, this allows the field at 
the center of mass of a system st tine t 0 to be expressed by superposing centrlhhtIons 
from each dmj at time r^/e earlier. 

It can be shown that the fiet effect of the retardation conditions on the source 
moving with uniform velocity is to make the total force appear to eome from the present 
rather than the retarded position or the source.When additionally the source ac¬ 
celerates, residual tangential components can arise in the force field and cause a 
degradation of the energy of an interacting sysfeat of masses, he examine this sit- 
fcatlstt in detail in the next sections. 

ll.C Held retardation. 

In I’art 1 we have seon how acceleration ®T a mass according to the Blrkhoff Theory 
modifies its static gravitational field, hhen a force acts on a mass causing it to 


Tidal frictions and the like are not considered in the present study. 


r 
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accelerate, It also acts to accelerate the negative equivalent mass of the gravita¬ 
tional field energy. The result Is summed up by equation (5), which holds down to 
the 1/c 9 order of residual gravitational effects. 

Craef has shown,as pointed out in the last section,that hjj must actually be 
treated as a retarded tensor potential. As mentioned, the net effect on the New¬ 
tonian potential h 0 Is to make the gradient «>f ti 0 at the origin appear to come from 
the present position of the sourde* to a very good approximation. The modified tensor 
potential however also Involves derivatives of the retarded velocity components of m. 
When expanded out in a Taylor series about t fee point t 0 - r/c, r being the distance 
to the center of mass of the syste», 

= -f*/0*v dt ♦(«•%! c 9 )d*M 0 /*i 9 .. m 

where Is the present vector velocity* We restrict oursolvcs to physical systems 
where v possewses continuous, well behaved doriv»tives|«»i % <e* Wo shall consider 
only terms down to the first order in !/«?• IV lWkt term of equation (61 how becomes! 

\/m * -l|h 0 « 4 /e* ♦ 2/f?Y* - fft^r/c Vv 0 A>t* (II* * 

The rate ‘af eftfcrgy dissipation true* I he dnnt tern in Itrul 

* ' fg = ♦ay 2 f?Y 2 - l)(mh 0 r/« 3 )% 0 *d e % 0 /dt f <m 

The minus' sign indicates the system dissipates power when *mh 0 “ -yMm/|U M ran he 
any mass at distance (at test fejative to the origin) wabjch interacts with is# 

11.1J Radial loir reaction In varions physical systems of interest. 

The result lends itself to an interesting pljysical interpretation: The gravita¬ 
tional force between m and M would be almost central except that the equivalent mass 
of their mutual field energy^which Is localized about nijlags m in space because of 
the above retardation condition on the velocity components. The direct gravitational 
force between that (negative) equivalent mass of the field energy and M then tries to 
pull m “backwards" and opposes tiie motion. 
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In the Ideal case of a mass m rotating lji a circle with h Q a a constant and with 
v 2 an invariant of the motion, the velocity Is perpendicular to the acceleration and: 

0 

v • F ■ -2yVr S - 1) (yjta/r) <ru/c) V (13) . 

• 8 -* 

Tljls is proportional to the square of the absolute acceleration and apparently has 
no angular dependencei Since the time average of the square of the acceleration Is 
tlte square of the acceleration, the result for a ring of mass orbiting around a cen¬ 
tral mass M becomes 

<V*M - -2Y 2 (SY 2 - 1) fyMr* /r)4rw/c) 3 w (14) 

S ring - ring 

whereas if the ring of mass holds itself together gravitationally and spins at the 
same time around its center of mass, 

v*f = -2 y 2 (3y 2 - l)iYm 2 (ruV/c 9 (15) 

lor art equal) balancing mass pole rotating in opposition to m, both v' and v are 
reversed and the power losses of the ««ntoin(jd masses are additive. 

Equation |13^ gives the energy dissipated in a unit time from an M-m planetary 
aystem» Evidently, the energy*is lost continuously and at a constant rate. The 
ratio of the power dissipated to the total negative energy ©f the M-m sf%t«w $s posi¬ 
tive so that the Ikinetie energy continues to increase and the planet spirals in 
towardfe Mi lor the Earth-Sun system, letting T be the present period of 1 year, 

dT/dt * -2(r(o/c) 3 toT (16) 

Since all the Values are known, immediate substitution and integration of the effect 
over 10 years results in the prediction that the year will be 10 hours shorter than 
how dUe to the gravitational reaction derived above. The present long term stability 
of atomic clocks falls Just slightly short of the- accuracy required for the measure¬ 
ment of vich an effect. Hotli m and M should radiate a power due to their own abso¬ 
lute accelerations* If r and a are the distances of m and M from the center of mass 







of tiie m-M system, mr = Ma; F m • v n * is found to equal * v^: and the sum equals: 


-4y (2y - 1)ym (ru ) /c 




The power dissipated by the mutual Interaction of m and M Is therefore M/m times 
the power either mass dissipates Individually.’ 1 ' As the ratio of M/m approaches i 
as In the case of the rotating dumbbell, the power dissipated in the mutual Inter¬ 
action approaches that dissipated by the masses individually* 

To calculate the self Interaction of an accelerating mass, let the center of 
mass of m be at Q(r,t) while P(r*, t Q - O' - r/c) is a non-accelerated point of an 
inertial space •Instantaneously" occupied by m. The field at P originating at Q 
is retarded t»|r a time O’* - Considering the Interaction of dm* with the static 

field of Mr we find the approximate |>ower dissipated in the self-interaction by 
integrating dm' mrer the entire mass distribution. This gives a net resistance 
to m*a (linear) autelepation ofl 


f s -Y 2 fBY 2 - llY« 2 (d 2 f 0 /dt 2 ll/c 3 (18) 


The power dissipated by the m-m dumbbell will be compatible with this result if 
we assume coherence will make the final result N 2 = 4 times the power from one of 
the masses.. 


See equation (13). 
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III. THE FAR FIELD 

III.A How does acceleratlon%ffect the energy density of the far gravitational field? 

In accordance with the findings of Tart I, only components^of gravitational field 
in the direction of absolute acceleration of source mass m will be relevant to first 
order calculation of the far field gravitational radiation intensity. This appears 
reasonable on grounds £hat field energy density has equivalent mass which reacts in 
line with the source's apparent acceleration as shown in equation (5). ® 

The Lagrangian function Moshlnski chooses (Postulate B.9) for the energy density 
of a Hirkhoff gravitational field reduces in free space in absence of other fields to: 

L = (-1/8) 7iYg 1J (dti^/a* 1 ) (19) 

r 

In this section we concern ourselves with the way the first order radial component 
**lr kjj changes along r and contributes to outward flow of gravitational radiation 
from the accelerating mass system. Thus we write 

L * (-l/8)nYg ij Oh lr /0x i Hdh lr /dJt J ) = (-l/4)TTY(Vh lr ) 2 (20) 

lr o 

where h^ B -h * -2 y h^/c* and we use this to get the self-interaction of the 
far field. 

The energy density in the far field t>y symmetry must propagate outward in the ra¬ 
dial direction with finite velocity c (Postulate B.5). In time dt, an amount of energy 
dE * L dV = Lc dt dA then flows out through the differential area dA. The total time 
rate of energy change within the volume enclosing the accelerating mass system is: 

dE/dt = j> A Lc • dA (21) 

Since dA = r 2 dfi where f) is the solid angle subtended by dA at distance from the origin 
r, only terms in the field energy density (20) whjcn go as 1/r 2 will have no singu¬ 
larities at infinity and will indicate the net gravitational radiation flow dependent 
only oh source parameters. We therefore exclude terms of higher order than 1/r 2 . 
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This exclusion will be Indicated by tne new symbol ■ which signifies that only 1/r 
dependent terms of the field and 1/r 2 dependent terms of the field energy density L 
are to be considered. 

111,11 Linearly and circularly accelerating mass. (Diagrams 1-2) 

© 

Let m, moving from rest at 0(0,0,0,0) with acceleration a along the x-axis, be 
at a point Q(t.,x,0,0) while P(t' ,x' ,y' ,z f ) is a non-aceeierated point of "absolute" 
space at which the field energy density is to be found. From classical mechanics, 

£mv 2 = max holds between the velocity, acceleration, and x coordinate of m at any 
time t. Thus W = ff/v. Since the y and z components of velocity are zero, (0) gives 

hj^ = -2Y 42 h 0 fv x /c)cos (r,v) (22) 

where eos(r,v) is the angle between OP and QP. The gradient of h is 

Vh Jr 5 (-2y 2 n 0 /c) cos(r,v) a/v (23) 

and equation (!!0) immediately gives L. Using this in conjunction with equation (21), 
we get finally* 

dt/dt g -| A (dfi/47i) (4y 4 ym 2 /c) (a/v) 2 cos 2 (r,v) = M/ay^yn^/c) (a/v) 2 (24) 

where if (ft is the angle between the projection of r on the x-y plane and the x-axis^ 
and 6 the angle between the r and tne z-axis, dflcos 2 (r,v) = sin 3 # cos 2 <t> d# d0. 

because overall momentum must be conserved for physical systems, the picture of 
an accelerating monopole of mass is too ideal. In reality, in3 1 = M3 2 where M moves 
In the opposite direction to counterbalance the motion of m. We now take m equal 
to M and find how this modifies the total far field potential. 

The velocities of the two mass points m and m are now equal and opposite, and 
h\ Postulate d.7 the total potential at a far field point P is: 

= -2Y 2 (h 0 /e) - Vj [ 2) )ros(r, V ) (25) 


ir 
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So far retardation conditions on the field interactions have been neglected, in 
which case the above potential would vanish. A little thought however shows the po¬ 
tential of the nearer mass to be slightly out of phase with that from the further 
mass due to the finite field propagation velocity c. The difference in phase of the 
two masses as seen at P is Just 2(x/c)cos(r,v). A previous relation (10) then gives 
to the first order I 


v< 2 > * v^ - 2(x/c)cos (r,v) (20) 

The total potential therefore does not vanish, and its gradient 

o e 

VHir % -4y 2 (h 0 3/c 2 )cos 2 (r,v) 

It follows immediately that: 

dE/dt § Lc • dA = f-4/ifr 4 Y(2m) V/c 3 (28) 

An equivalent line of reasoning gives similar results to the above for the case 
of circularly accelerating mass. Again, calculation of h lp requires that the velocity 
of the circling mass be projected onto the radius vector to the far field point P. 

The angular dependence of the h^ p field becomes sin0 sin(<£ - wt) where u is the an¬ 
gular velocity of rotation and 6 and <£ retain the same meaning as above, lor the 
case of the rotating dumbbell, retardation conditions are accounted for in the phase 
difference in the sin(<£ - wt) term. The total potential at P then becomes: 

• H = -4y 2 (h n x*3/c 2 ) sin 2 # cos 2 ($ - tot) (29) 

Ir 0 

The total power dissipated can then be gotten as above by computing the gradient of 
H^ p and from it L, then integrating L over the appropriate enclosing area. This 
of course gives back exactly the result in equation (28). 

There is another point of even more interest in the expression for in this 
case of circularly accelerating mass. Namely if represents contributions from 
opposite mass elements of a rotating ring of mass, integration over all the mass in¬ 


is approximately: 

(27) 





volves the term oos*(^: - u>t)dfi integrated from 0 to it. But this integral dqes. not 
® vanish, and therefore Hj r at a far field point although becoming time independent 
in this case does not vanish. Neither does its gradient or the outward flowing 
gravitational field energy density L! Integration of L over the ; enclosed area con¬ 
taining the rotating mass rin g then gives a steady flow of energy out of the system: 

dh/dt g f-8/15) y 4 yM V/c 3 (30) 

Here M is the total mass of the ring. 

Equations (28) and (30) should be compared with (17) and (15) for the radiation 
reaction for a mass dumbbell and ring derived in Part II. The agreement is fairly 
close, but not exact numerically. The shortcoming may derive from not considering 
other components of h^j in the second*order computation involving the retardation 
conditions, or it may be irremovable in the same way that the electromagnetic mass 
of the electron does not account for its entire mass. 

There is an interesting physical interpretation of what happens when m acceler¬ 
ates. In the Birkhoff Theory acceleration of a source mass has been seen to modify * 
its gravitational field. One can imagine the field acquires a new energy density 
which propagates out from the mass at c. One might look on the induction field action 
also as propagating out at c in the static field case. However the net energy 
"flawing out M of such a E field by (20) goes as 1/r and thus rapidly drops to zero 
with increasing r. The net energy flow from the residual accelerated part of the field 
has ho r dependence. It therefore gives a net energy loss at distances far from the 
system. 

This type of gravitational radiation cannot be thought o^ in terms of oscillating 
multipole fields. It appears solely a result of the physical acceleration of an 
Interacting mass system. 

III.C Coherence. 

Hhen N masses make up each pole of a spinning dumbbell, (28) tells us the power 
.radiated is X £ times the single m-m dipole. This is the case- for radiation from N 
coherent sources. 17 lor most macroscopic^physically realistic systems, the wavelength 
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X e 3nc/b) is much greater than the physical dimensions. Ihe coherence is a 

system 

result of residual radiation components emitted in phase with each other. 

Next, one might naturally ask what is the gravitational radiation loss for a 
Ntn-m dipole spinning about its center of mass. In this case m undergoes an acceler¬ 
ation a 0 and Nm an acceleration a Q l/N about the center of mass. (28) shows each 
pole radiates proportional to its gravitational field energy and the square of its 
cwn acceleration. The two masses taken with their respective accelerations are thus 
found to radiate equally. Due to coherence, the net power is ^Just 4 times that 
expected from the single accelerating pole alone. 

So far the mutual gravitational field energy of m and Nm has not been taken into 
account. Assuming the mutual field energy is localized about m for reasons explained 
in Parts I and II, the accelerational interaction for large N dissipates a power 
equal to N times that dissipated by the single accelerating pole alone. In general 
an "amplitude'' can be defined proportional to the square root of the gravitational 
field energy and to the absolute value of the total acceleration. Because of coher¬ 
ence the power is the square of the total amplitude, but for large N the cross terms 
and the separate contributions from m and Nm can be neglected with respect to the 
radiatioh term caused by the mutual gravitational energy. Equivalently, the single 
pole m can by thought of as interacting separately with each one of the N masses of 
the opposite pole, each pair contributing incoherently to a net power again approx¬ 
imately N times that from the single accelerating pole alone. 

An m-M gravitational system can therefore lose energy into its far field at a 
rate proportional to the mutual gravitational potential energy of the system. Such 
an-interpretation is consistent with the estimates of radiation reaction found in 
Part II. The radiation reaction dissipation (13) has therefore been correlated with 
ttie far field gravitational radiation in the framework of the Birkiioff Theory. As 
mentioned previously, as m approaches the size of M, i.e. m = M, the two interpre¬ 
tations definitely become equivalent, and the total power radiated goes as the mutual 
potential interaction energy of the masses. 







IV. TWO PICTURES 01 GRAVITATIONAL RADIATION LOSS 


IV.A Integration of the preceding views of radiation from accelerating masses. 


It may prove helpful to integrate the particular examples of tlie preceding sec¬ 
tions into a general view of gravitational radiation in flat spacetime. The perfect 
fluid* of BirkhofT, though purely an auxiliary concept and not essential to the Birk- 
hoff Theory, offers a formal connection between the gravitational and inertial forces 
and will be used nere to aid in linking the complementary pictures of direct inter- 
particle radiation reaction forces and the far gravitational radiation field intensity 
The equation of state of the Uirkhoff perfect fluid, whose internal mass-energy 

and momentum is characterized by the tensor T 1 ^ = T^ 1 and gravitational field by 
1 1 2 

tg (r,ct) at each point, is defined bj p = p 00 c where the pressure and density are 
scalar world invariants. 

T*'J describes solely the state of the "true 1 * mass. It therefore vanishes every¬ 
where outside the world tube of a mass particle m. t^ (r,ct) describes all comple¬ 
mentary gravitational field effects. With only gravitational forces acting, the 

divergence of the sum of and t*J(r,ct) must vanish to ensure equilibrium of 
e ® 

forces* When the proper test-energy density as here is assumed a priori constant 
along the world line of m, th® total energy density generally can fluctuate between 
kinetic and gravitational. In absence of gravitational fields, this provides the 
standard continuity relations for the perfect fluid, while in absence of "true" 
masses we obtain the conservative character of the static gravitational field and 
continuity of any radiation fields which may occur (although equivalent mass in the 
field may also be a source of radiation). 

Integrating the divergence of T 1 ^ + t^(r,et) over space, and applying the theo¬ 
rem of Gauss in 4 dimensions, <(, [T 1<X ],dA = 0 since 1’ 10£ = 0 outside the world tube and 
J A A 

J v Ydt“ dV 0 = -J v y^ t“(r,ct)dV - j {ctg a (r,ct)} A dA a (31) 


By direct integration, the left term is found to be the time component, of inertial 
Minkowski 4-force (yd/rit)ymc Z , the time rate of change of energy in the mass system. 
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IV.B® Radiation Reaction Picture. 

* 

The time rate of change of total inertial energy can be seen to be equal to the 

time component of llirkhoff gravitational 4-force multiplied by c. The sum of the 

i Ik 

last, two terms of (31) is then by definition just me B ^ u J u . This represents the 
rate at which external gravitational forces acting on m do work. Because of the con¬ 
servative character of the gravitational induction field, this therefore expresses 

the residual gravitational radiation reaction when averaged over a full cycle or 
ilk 

long time interval. Note that me B ^ u J u depends only on relative parameters, i.e. 
mutual energy, relative velocity and acceleration, etc., all in a framework of flat 
spacetime. Thus we have a direct interparticle calculation of radiation damping or 
friction as opposed to the complementary picture which handles gravitational radiation 
solely on a basis of field interactions. 

> <* 

IV.C Direct Held Interaction Picture. 

^Besides calculating directly the self or external field reaction of each mass 
element dm and summing over all elements to obtain the total time rate of change of 
energy, a far field radiation intensity, which should yield a compatible result, can 
be obtained by a complementary method. Let us say the energy of a gravitational ra¬ 
diation field at distance r from m propagates in the direction of r with fundamental 
velocity c. The energy crossing an area dA in time dt will be found at far distances 
in a volume c dA dt. Using Moshinskl’s field energy density,Lagrangian function, 
the total power dissipated over an enclosing surface will be: * 

dL/dt ” Lc • dA = {(c/8ny)g 1 ^(9h ,nn /dx 1 ) Qh^/Qx^ } • dA (32) 

Tills result can be compared with the selection of the Lagrangian energy density for 
the electromagnetic field and with the corresponding form of the electromagnetic 
Poyntlng vector. 4 
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The first term on the right in relational) vanishes for our purposes because, 

a) the time average of the conservative induction field is independent of time, and 

b) the local energy density of the gravitational radiation field itself can be 
assumed in tne case of "steady*' radiation processes to be constant on the time av¬ 
erage. The last term describes the total flow of gravitational radiation energy 
through a bounding surface and will be taken to be independent of the distance r 
from the center of tne radiating mass system. Under these conditions, the quantity 
in brackets, {ct^ a (r,ct)}, compares with the quantity in brackets in the previous 
relation (32), and both of these take the form of a gravitational Poynting-type- 
vector describing the radiation flow, a function in general of the distance r and 
the polar angles 8 and <£ 0 - wt. 

The dual aspect of our radiation picture is now apparent. In the case of the 
far field, the angular and 1/r 2 dependence of the gravitational radiation inten¬ 
sity will enter through relation (32), whereas in the direct interparticle inter¬ 
action of mass m with self or external gravitational forces, no necessity appears 
for evoking the angular dependence of the emitted radiation. Resolution of this 
possibly significant difference in our radiation pictures would require more scru¬ 
tiny than applied here, and remains a problem of certain interest. A point in fa¬ 
vor of these pictures of gravitationaleradiation is that both otherwise seem com¬ 
pletely consistent with the Birkhoff Theory, which indeed was originally intended 
- 

to describe completely all gravitational effects. One must only be aware that 
more refinements, such as the retardation conditions imposed by Graef, 10 are required 
in the treatment of effects of higher than 1/c 2 order. 
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V. REACTION OF CENTRIFUGAL FORCES ON FREELY SPINNING MAS8ES. 

V.A Interaction of a freely spinning mass with the universe. 

Interaction of the gravitational field of an accelerating mass element m In a 
distribution of elements whose center of mass remains at rest with respect to a large 
remote external mass M can be calculated directly from equation (4). Noting the 
only non-vanlshlng component of v^ In this case Is v 1 » 1, and If F g Is the relativ¬ 
istic 3-vector force which m exerts on M at a distance 1( so great that 1/n* and 
higher order terms can be neglected, the last term of (4) becomes: 

F 5 2Y 4 [1 ♦ (v/c)*] ( Y lfe«/Il)8*/c’ (33) 

ft - 

Here v denotes the velocity of the element m relative to Its system's center of mass 
(and hence In this case to M). 

It may be well before going on to indicate the physical meaning of (33). 

Previous arguments of 0^0^ suggest that In the Hlrkhoff Theory the limiting value 
of YMm/hfas for example for the universe)Is £mc 2 . Were this Indeed the case, 

F g g y“[ 1 ♦ (v/c)*]m3+ • (34) 

that Is, ttiere Is a gravitational Interaction with the universe approximately equal 
to the Inertial (centrifugal) lorca on m. Indeed, this expresses Mach's contention 
that the so-called Inertial force Is actually of gravitational origin. The mutual 
gravitational field energy of m with respect to M again appears to be localized 
about m and to undergo the same acceleration. 

Going on, we now ask what residual effects may occur In the Interaction of the 

* 

mass elements of a freely spinning mass with the universe. Because the acceleration 
of an element m must be referred to the center of spin, by equation (10) it must be 
treated as a retarded quantity. This results In the Taylor Series expansion analo¬ 
gous to equation (11) wnlch gives a residual component of Interaction force In line 
with the vector velocity of m. Neglecting all force reaction effects on the rest 
energy me of m, which has been assumed constant along the world line, and using 


'Orbit Theory", Proc. bymp. App. Math. Vol.IX, p.171, (1957) 
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the results of equations (0), (11), and (33), the kinetic energy of each mass element 
Is found to decrease at a rates 

dK.L./dt = *v = -4y 4 (K.E./c 2 )(r/c)v • d*v/dt 8 (86) 

where r Is the spin radius of m and K? E. its kinetic energy. Since the product 
of v and its second time derivative is the same sign and a colinear product for an 
isolated, freely spinning system, the losses from each mass element are apparently 
additive and contribute to a net energy decay from the entire system. The "freely" 
spinning mass therefore will slow down due to its interaction with all masses in 
the universe. 

This is strictly a gravitational type interaction. Physically, the net power 
dissipation in this case appears to be non-oscillatory and independent of direction. 
If the universe by definition cannot itself spin, the entire angular momentum re¬ 
action force must set back on tne freely spinning mass system causing a real slowing 
down. 

* 

V.B Self interaction of mass elements in a freely spinning system. * 

Equation (34) suggests that centrifugal forces can be interpreted in the Birk- 
hoff Theory in terms of direct gravitational interaction with the universe. These 
forces according to postulate 11.5 propagate^at finite velocity c. One can therefore 
ask what residual effects arise in interaction of oppositely paired elements in a 
freely spinning mass distribution. 

Consider Diagram 4. We ask whether the line of force between oppositely paired 
mass elements passes througn the center of mass. If it did, it could be assumed no 
tangential components of force were acting on either mass. However we shall show how 
acceleration enters the picture to preclude the trivial result. 

It takes a time 2a/c in the proper rest frame of the spinning system for a force 
disturbance from one mass element to propagate to the other. This has led or perhaps 
misled some authors*' 1 ^® to assume the existence of 1/c order tangential forces within 
the system. However it is not hard to see that .in the proper frame the oppositely 
paired elements are at rest with respect to each other and no such components actu¬ 
ally arise. However the propagation of centrifugal forces is invariant only with 
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respect to velocity In a flat spacetime, not to acceleration or its rates of change, 
as evidenced by tlie bending of light in a gravitational field. With respect to the 
propagating centrifugal forces which are not acceleration invariant,the opposite mass 
element in time 2a/c accelerates a distance in the tangential 6 direction of l/fla(2a/c) 8 . 
The angular discrepancy which must appear in the centrifugal force is therefore 
4/3(aa 2 /c 3 ) so that a 1/c 3 order component of that force does lie in the 6 direction. 

When this component is dotted with the velocity and the relativistic factors taken 
into account by transforming the process out of the proper frame of the spinning 
mass, the result for circular motion is Just 2/3 Of the result found in equation (36). 
This suggests the two treatments, V.A and V.B, are very nearly equivalent. 

Physically, in the rest frame of the spinning system the centrifugal forces would 
appear to diminish with time. This describes a new process which can be called the 
centrifugal damping of spin motion. m 

lor the spinning point dipole of radius 'a', equation (35) gives: 

dE/dt = 2F^‘ v = -4Y 4 mv 2 (aw/c) 3 w (36) 

and the proportional decay of kinetic energy for low velocities is: 

d(K.E.)/dt - + 2d(i)/dt * -4 (aw/o) 3 w (37) 

KX a) 

It is particularly interesting that spontaneous loss of rotational kinetic energy, 
similar to the losses predicted from the gravitational systems in Parts'll and III, 
does not appear to cancel when integrated over a mass such as a rotating sphere or 
cylinder whose external multipole moments do not vary with time. Since the "grav- 
itational radiation" reaction in the case of the universal interaction is propor¬ 
tional directly to m, tiie contributions of £ach mass element to the total power 
dissipated are additive and equation (36) can be generalized to: 

( dK/,lt > general radiation reaction = “{» &Y 4 (r^/c) % V dm (38) 

This expression sets a lower bound on the rate "freely" spinning mass systems slow down. 
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A lucid interpretation attaches to the separate situations of purely gravita¬ 
tional and centrifugal radiation damping. Physically, in (13) picture v/c compon¬ 
ents of a negative outward-flowing "g-field" fluid acting against the negative 
equivalent mass of the m-M field. In (30) picture v/c components of positive out- 
ward-flowing "centrifugal field" riuid to act against the positive equivalent mass 
of the kinetic energy. A little reflection will show both effects go in the same 
direction, causing a net rate of energy loss from the accelerating mass system in 
agreement with the minus signs which appear in these formulae. 

Gravitational or/and centrifugal damping of a freely spinning mass system's 
motion therefore results in a positive outward flow of 1/c order gravitational 
radiation energy,accompanied by a genuine decrease in the total energy of the source 
system, i.e. a genuine decay of the motion. This is in contradistinction to the 
net erfect of l/c“ order gravitational quadrupole radiations predicted by Einstein 
and others, 3 whereby a negative energy mysteriously appears to be radiated outward 
while the total internal energy of the mass system increases. (See, for example, 

A. Peres, II Nuovo Clmcnto , XI, 1 March 1059, p. 053). 

V,L Expe rimental magnitude of the predicted 1/c 3 order gravitational radiation. 

lor physically realistic, laboratory-sized, symmetrically rotating masses, 
acceptance of equation (38) leads to a proportional decay of angular frequency for 
freely spinning masses of magnitude: 

* (-2Y 4 w 4 /r 3 ) j m r B dm/J ni r'dm * -k(a/c) 3 w* (39) 

1* depends on the geometrical moments of the mass, being 8/7 for a right cylinder 
rotatit* about its axis, 1 for a thin rod rotating about an axis through its center 
normal to its length. The minus sign indicates that the spin frequency decreases. 

Hith a view to experimental feasibility, for a steel cylinder spinning near 
the limits of tensile strength at velocity au> * 10 B cm/sec, the frequency decay 
ratio depends on u\ Let a * 10 3 cm with u> on the order of lo" radians/sec. 
dw/dt is then found to be on the order of * radians/sec/sec. The effect is thus 
on the order of ft parts in in 9 on frequency. 
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It. seems within possibility slmuljtaneouslj to generate and measure such quanti¬ 
ties using a type of magnetic suspension and rotation apparatus developed by Beams. 

By suspending the rotating mass in ultra-nigh vacuum and making use of suspension 
symmetry, virtually all stray frictions are impressively reduced. Moreover, magnetic 
and other residual drags appear to vary proportional to much lower powers of the 
frequency. 20u 

The angular frequency and its rates of change are capable of being detected 
with great precision by comparison with a good frequency standard. The experimental 
techniques appear hard but straightforward. 29 * 3 On the other hand, tests for the 
gravitational quadrupole radiation of Einstein 3 , even if reinterpreted in the sense 
of sections V.A and V.B, are of an (aw/c) 2 lower magnitude and far beyond the range 
of present instrumentation. 
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